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ABSTRACT 

As part of an ongoing study to determine the distribution of orbital parame- 
ters for massive binaries in the Cygnus 0B2 Association, we present the orbital 
solutions for two new single-lined spectroscopic binaries, MT059 (08V) & MT258 
(08V), and one double-lined eclipsing binary (Schulte 3). We also constrain 
the orbital elements of three additional double-lined systems (MT252, MT720, 
MT771). Periods for all systems range from 1.5-19 days and eccentricities range 
from 0-0.11. The six new OB binary systems bring the total number of multiple 
systems within the core region of Cyg OB2 to 11. The current sample does not 
show evidence for a "twin-heavy" binary distribution. 

Subject headings: techniques: radial velocities — (stars:) binaries: general - 
(stars:) binaries: spectroscopic — (stars:) binaries: (including multiple) close - 
stars: early-type — stars: kinematics — surveys 



Introduction 



The formation of massive s tars is the subject of extensive debate (jKrumhola 12005 



Bonnell et al.lll998l : lLarsonll200ll ). Current formation mode ls are based, in part, on the ac- 
cumul a ted statistic s of massive star bina r y syst ems, such as Evans et al.l (120061 ). iMermilliod 
(119951 ). iGied (119871 ). and iGarmany et al.l (jl980l ). The extent to which these statistics can be 
generalized is limited, however, because some studies only sample small numbers of binary 
systems while others include large samples of systems from radically different environments 
wi th varying ages. A re c ent ad aptive optics imaging study and modeling analysis of Sco OB2 
by iKouwenhoven et al.l (120071 ) has attempted to circumvent the selection effects of previous 
smaller studies in order to obtain the true binary fractio n of this cluster. They find a binary 
fraction for B and A-type stars between 70% an d 100%. iKobulnickv &: Fryer! (120071 ) utilized 
Monte Carlo simulations applied to the data of iKiminki et al.l (120071 . hereafter Paper I) to 



1 Dept. of Physics & Astronomy, University of Wyoming, Laramie, WY 82070 
2 Dept. of Physics, Lehigh University, Bethlehem, PA 18015 



-2- 



find a true binary fraction of / > 70% for O and B stars in Cyg OB2. Unfortunately, the 
true binary fraction of massive stars in young clusters is poorly known, but it is one of 
several preserved parameters providing information about the formation of stellar systems. 
Additional parameters preserved from the formation of young binary/multiple systems are 
the eccen tricity, separ ation, period, and angular momentum (in the form of rotational ve- 



locities) (jLarsonll200ll ). Understanding the initial conditions that lead to the formation of 
massive stars is aided by reliable statistics on a sample of binary systems sharing a common 
formation history. As a step toward understanding the characteristics of massive binaries in 
Cygnus OB2, we present a study of six new massive systems from the radial velocity survey 
of 150 OB stars begun in Paper I. The OB stars are probable memb ers of the Association 



and li e within ~15 pc (at d = 1.7 kpc) of the cluster center surveyed by lMassey &: Thompson 



(Il99ll ). The long term goal of this study is to provide an accurate distribution of binary 



orbital parameters for massive stars with an assumed shared formation history. 

§ 2 provides observational details of the new spectroscopic datasets. § 3 discusses the 
measurement of radial velocities, the search for periods in the radial velocity data, and the 
determination of orbital elements via radial velocity curve fitting. §4 discusses the orbital 
solutions to the single-line d binaries (SBls), MT059 & MT258 (notation in the format of 



Massey &: Thompson! 119911 ). and the double-lined binary (SB2), Schulte 3 (Cyg OB2 no. 3). 
§ 5 constrains the orbital elements of three additional SB2s, MT252 & MT720, and MT771. 
Finally, § 6 summarizes the results of the survey to date, including the total number of OB 
binaries uncovered in the Cyg OB2 core region. 



2. Observations 

Paper I details the observations of this survey through 2005 October. We have obtained 
seven additional datasets with the Longslit spectrograph on the Wyoming Infrared Obser- 
vatory (WIRO) 2.3 m telescope and the Hydra spectrograph on the WIYNl^J 3.5 m telescope. 
Tabled] lists the observing runs at each facility and the corresponding spectral coverages and 
mean spectral resolutions. Instrumental setups at each facility for the first six runs were 
consistent with previous observations described in Paper I. 

Observations at WIYN took place over seven nights on 2006 September 8-11 and 2007 
July 3,4, & 5. We used the Hydra spectrograph with the Red camera, 2" blue fibers, and 
the 1200 1 mm -1 grating in second order to obtain three 1200 s exposures in each of three 
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fiber configurations (1500-1800 s on the second run) yielding a maximum signal-to-noise 
ratio (S/N) of 80:1 (2006) and 130:1 (2007) for the brightest stars. The spectral coverage 
was 3800-4500 A at a mean resolution of R ~4500. Copper-Argon lamps were used between 
each exposure to calibrate the spectra to an RMS of 0.03 A (2 km s" 1 at 4500 A), and 
the typical resolution was 1.0 A FWHM at 3900 A and 0.82 A FWHM at 4400 A. Spectra 
were Doppler corrected to the heliocentric frame and checked against the radial velocity 
standards HD182572 (G8IV), HD187691 (F8V), and the minimally varying, amply sampled 
star MT083 before comparison to previous datasets. Both runs were plagued by intermittent 
clouds. 

Observations using the WIRO-Longslit spectrograph with the 600 1 mm -1 grating in 
second order took place over 12 nights from five runs of varying length (2006 June 16-20, 
2006 July 15-16 & 19, 2006 October 7, and 2007 June 28-30). Exposure times varied from 
600 s to 1800 s depending on weather conditions to obtain a maximum S/N of 150:1 for 
the brightest stars. The spectral coverage was 3900-5900 A. Copper-Argon lamp exposures 
were taken after each star exposure to wavelength calibrate the spectra to an rms of 0.03 A 
(1.9 km s _1 at 4800 A), and the typical spectral resolution was 2.9 A FWHM at 4100, 2.4 A 
FWHM at 4800, and 2.8 A FWHM at 5700 (2006 October and 2007 June had a more 
uniform resolution of ~2.5 A FWHM across the chip). Spectra were Doppler corrected to 
the heliocentric frame and checked against the radial velocity standard HD187691 and the 
minimally varying, amply sampled star MT083 before comparison to previous datasets. The 
2006 June and July data were affected by intermittent clouds and an abnormally high dark 
current arising from a camera cooling problem. 

We also obtained 8 nights of spectra using the WIRO-Longslit spectrograph with the 
1800 1 mm -1 grating in first order (2007 August 28-September 4) to examine the Ha and He I 
absorption lines in suspected SB2s. Exposure times varied from 600 s to 1800 s depending 
on weather conditions and yielded a maximum S/N of 200:1 for the brightest stars. The 
spectral coverage was 5550-6850 A and Copper-Argon lamp exposures were taken after 
each star exposure to wavelength calibrate the spectra to an rms of 0.03 A (1.4 km s _1 at 
6400 A). The typical spectral resolution was 1.5 A FWHM across the chip. Spectra were 
Doppler corrected to the heliocentric frame and checked against the radial velocity standard 
HD 187691 and the minimally varying, amply sampled star MT083 before comparison to 
previous datasets. All new datasets were reduced using standard IRAF reduction routines 
as outlined in Paper I. 
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3. Data Analysis and Orbital Solutions 



We obtained radial velocities, V r , for the new observations via the IRAF cross- correlation 



task XCSAO in the RV SAO package (IKurtz et al.lll99ll ). using a model stellar atmosphere 
(ILanz fc Hubenyl 120031 . TLUSTY) of the appropriate ef fective temperat ure and gravity as 
discussed in Paper I. We used the method outlined in iMcSwainl (120031 ) to determine the 
orbital parameters of each system. We obtained an estimate of the orbital period through an 
IDli program written by A. W. Fullerto n which makes u se of the discrete Fourier transform 
and CLEAN deconvolution algorithm of iRoberts et al.l (119871 ). The strongest peaks in the 
power spectrum of each star were examined by folding the data at the corresponding period 
and inspecting the V r curve visually. Orbital elements were then procured by using the 
best period as an init i al est imate in the nonlinear, least-squares curve fitting program of 
Morbey &: Brosterhusl (119741 ) . The best solutions were attained by manually varying the 
initial guesses of key orbital parameters until a minimum in the standard deviation (rms) 
was found. Weights for each point were assigned as the inverse of the la cross-correlation 
velocity error, where error within the XCSAO task is calculated using, 



3w 



1+r 



(1) 



In this equation w is the FWHM of the correlation pe ak and r is the rat io of the correlation 
peak height to the amplitude of antisymmetric noise (IKurtz et al.lll99ll ). 



We also included an additional 10 km s _1 error added in quadrature to the 1 a cross- 
correlation errors for all spectra from the 2006 June and 2006 July datasets. A varying offset 
of up to ~10 km s _1 was measured for the radial velocity standard HD187691 and the amply 
sampled, minimally varying system MT083 during these runs. We have not been able to 
ascertain the source of this velocity anomaly, and we did not encounter this problem during 
other runs. 

The complete list of orbital elements for each binary system appears in Table [2j Listed 
within the table are the period in days (P), eccentricity of the orbit (e), longitude of peri- 
astron in degrees (uj), systemic radial velocity (7), epoch of periastron (To), primary and 
secondary semi-amplitudes (K\ & K 2 ), adopted or calculated minimum primary and sec- 
ondary masses in solar masses (Mi & M 2 ), primary and secondary mass functions in solar 
masses (f(m)i & f(m)-i), spectral classifications from this survey (S. C.\ & S.C2), the mini- 
mum primary and secondary semi-major axes in solar radii (a\sin i & a2sin i), and finally, 
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the rms of the fits {rms\ & rms2) 



4. Orbital Solutions to MT059, MT258, and Schulte 3 



4.1. The SB1 MT059 



We have accumulated 45 observations of the 08V star, MT059 (Cyg 0B2 No. 1), includ- 
ing six at Lick, two at Keck§|, eight at WIYN, and 29 at WIRO. A total of 44 spectra were 
included in the period search after eliminating one low S/N Lick observation. The singular 
peak in the CLEANed power spectrum occurs at a period of P = 4.852 days. The power 
spectrum also contained a small amount of low power noise but no potential aliases. Figured] 
shows the V r curve and best-fitting orbital solution with a period of P = 4.8527 ± 0.0002 
days (the three most discrepant points in the figure belong to the 2006 June & July runs and 
were given less weigh t to improve the fit ). Given this period, an assumed primary mass of 
Mi = 21.4 ± 0.6 M jMartins et all 120051 ) . and a semi-amplitude of K x = 72.3 ± 2.3 km s -1 , 
the solution yields a minimum secondary mass of M 2 = 5.1 ± 0.3 M (where the mass error 
is derived from the error in the mass functio n). This corresponds to a mid B star or earlier 
(interpolated from iDrilling fe Landoltl I2OO0I ) and a mass ratio q = M 2 /M x > 0.24 ± 0.02. 
In the best spectra having S/N of 150:1 (WIRO-Longslit), we should be able to detect a 
secondary star's spectral features at velocity separations of (1 + l/q)K\ if the component 
luminosity ratio is larger than L 2 /Li ~ 0.2 (luminosity ratios have been chosen based on 
a set of spectra combinations at various luminosity ratios and line separations according 
to spectral type and theoretical mass ratios). The absence of spectral signatures from the 
secondary limits the luminosity ratio to L2/L1 < 0.2, which, in turn, implies a secondary 
later than B1V and an u pper limit on the mass of M 2 < 13.8 ± 0.6 M (interpolated from 
Drilling fc Landoltl |2000| ). Thus, the secondary mass is constrained to 5.1 < M 2 < 13.8 M 
and the mass ratio is 0.24 ± 0.02 < q < 0.64 ± 0.03. The upper limit on the secondary 
mass also requi res i > 26°, howev er this is a loose requirement given the uncertainty in 
O star masses (IMassey et al.ll2005l ). Additionally, the solution yields a minimum primary 
semi-major axis of aisin i = 6.9 ± 0.2 R , an eccentricity of e = 0.11 ± 0.04, a longitude 
of periastron of uj = 212 ± 17°, and a residual rms = 11.0 km s -1 . The residuals show no 
indication of periodicity and have a Gaussian distribution. 



The ephemerides for MT059 and MT258 are listed in Table [3] and include the star 
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designation, date (HJD — 2,400,000), phase (0), measured radial velocity (V r ), cross- 
correlation la error (la err), and the observed minus calculated velocity (O — C). 



4.2. The SB1 MT258 

We observed another 08V star, MT258, 42 times, including two at Keck, six at Lick, 
13 at WIYN, and 21 at WIRO. The six spectra obtained at Lick and one spectrum ob- 
tained at WIYN yielded low S/N and unreliable velocities with large uncertainties and 
were subsequently removed. The strongest signal in the CLEANed power spectrum for 
MT258 corresponds to a period of P = 14.66 days. In addition to some low power noise, 
the power spectrum also contains an additional relevant signal corresponding to a period 
of 0.49 days. Insufficient variations over the course of a night coupled with a visual in- 
spection of the folded V r curve reveals this period is an alias. The V r curve and best fit- 
ting orbital solution shown in Figure [2] (rms = 6.7 km s _1 ) corresponds to a period of 
P = 14.660 ± 0.002 days and an orbit with an eccentricity of e = 0.03 ± 0.05. Addition- 
ally, the orbital solution provides a semi-amplitude of K\ = 41.2 ± 1.7 km s -1 , a primary 
semi-major axis of aisin i = 11.9 ± 0.5 R , a longitude of periastron of uj = 68 ± 71°, 
and a minimum seco ndary mass of M 2 > 4.1 ± 0.2 M , assuming a 21.4 ± 0.6 M primary 



(IMartins et al.l 120051). Mo > 4 .1 ± 0.2 M corresponds to a late-to-mid B star (interpolated 



from iDrilling fc Landoltll2000i ) and a mass ratio q > 0.19 ± 0.01. Given the comparable S/N 
between MT258 and MT059 spectra, we should be able to detect a secondary star's spectral 
features at velocity separations of (1 + l/q)K 1 if the component luminosity ratio is larger 
than L2/L1 — 0.4. The absence of spectral signatures from the secondary and the semi- 
amplitude of K\ = 41.2 ± 1.7 km s _1 limits the luminosity ratio to L2/L1 < 0.4, which, in 
turn, implies a seconda ry later than BOV and an upper limit on the mass of M 2 < 17.5 ± 1.0 



M (interpolated from IDrilling &: Landoltll2000l ). Thus, the secondary mass is constrained 



to 4.1 ± 0.2 < M 2 < 17.5 ± 1.0 M and the mass ratio is 0.19 ± 0.01 < q < 0.84 ± 0.05. The 
upper limit on the secondary mass also loosely requires i > 18°. 



4.3. The SB2 Schulte 3 



Schulte 3 is lis ted as a probable binary system on the basis of Ha variability by 
(120021 ). and we report here the double-lined nature of this system for the 



Harries et al. 



first time. It is a pr obable member of Cyg OB2 , but it lies on the outskirts of the clus- 
ter core surveyed by iMassev fc Thompson! (Il99ll ). The spectral type of this star, 08.5IIIf 



(IMassey &: Thompson 



1991 



09: in the Simbad database) is not entirely certain, and there- 
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fore the adopted absolute magnitude from iMartins et al.l (120051 ) is also less certain. The 
spectrum does carry the same diffuse int erstellar band absorption (at 4428 A and 450 f A) 
characteristic of known Cyg OB2 stars (IHanson 1 120031 ) . Preliminary results of a concur- 
rent photometric survey of thi s cluster show this bina ry is also an eclipsing system with an 
estimated period of 4.75 days (jKinemuchi et al.ll2008l . in prep). 

We obtained 9 exposures of this star between 2000-2007 with Hydra at WIYN, one 
exposure on 2000 September 19 at Keck, and eight exposures between 2007 August 28 and 
September 4 at WIRO with the Longslit spectrograph. The five exposures obtained at WIYN 
from 2006 September 8-11 show a progression from deblended, double-lined profiles to a 
single, slightly asymmetric profile and back. Figure [3] shows this same deblended/blended 
line progression in the 2007 August /September WIRO data. Using 14 of the 18 observations 
obtained, we deblended the He I AA4471, 5876 A lines by fitting simultaneous Gaussian 
profiles with the SPLOT routine in IRAF to obtain rough V r measurements. We repeated this 
method 10 times while varying the baseline region used to define the continuum each time. 
We then used the mean of each component's list of Gaussian centers to compute the velocity 
and adopted the rms of this list as our one sigma uncertainty. This results in an average 
uncertainty of < 8 km s -1 for the secondary and ~ 20 km s _1 for the primary. We observe 
a maximum line separation of ~370 km s -1 for the He I lines. Emission in Ha associated 
with the primary makes it difficult to compare helium with hydrogen velocities, but the 
maximum measured V r separation for He I is consistent throughout the data. Based on a 
visual inspection of the spectra in Figure EH a period of 4.63 days from the CLEANed power 
spectrum, and the preliminary results of iKinemuchi et al.l (120081 . in prep), we estimated a 
period of 4.7 days. It should be noted however, that the CLEANed power spectrum also 
contains a stronger signal corresponding to a period of 1.29 days, a 1 — v alias (where v 
is the signal frequency in days^ 1 ). However, a period of 1.29 days does not produce a 
sinusoidal folded V r curve for the primary or agree with the other period estimates. We 
attribute the power of this signal to the considerably smaller number of data points used in 
the period search. The best orbital solution yields a period of P = 4.7464 ± 0.0002 days 
for the secondary. The period, eccentricity, epoch of periastron, and systemic velocity were 
applied as fixed parameters to obtain the solution for the primary. The residuals and errors 
are considerably higher for the primary and likely a result of emission in the helium line cores 
during a portion of the system's orbit. The V r curve and orbital solution for Schulte 3 is shown 
in Figure HI The filled symbols correspond to the primary V r measurements and the unfilled 
symbols correspond to the secondary V r measurements. With a low eccentricity of e = 
0.070±0.009 and semi- amplitudes of K 2 = 256.7±2.4 km s" 1 and K 1 = 113.2±14.5 km s -1 , 
the solutions yield a nearly circular orbit and a mass ratio of q = 0.44 ±0.08. The calculated 
semi-major axes are aisin i = 7.4 ± 0.9 R and o^sin i = 16.7 ± 0.2 R (i.e., a separation of 
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asin i = 24 ±0.9 R©), and the calculated lower limit masses are Misin 3 i = 17.2 ±0.3 M and 
M 2 sin 3 z = 7.6 ± 1.4 M . The ephemeris for Schulte 3 in Table H] lists the star designation, 
date (HJD— 2,400,000), phase ((f)), measured radial velocities (Vi & V 2 ), and the observed 
minus calculated velocities (0\ —C\ & O2 — C2). Errors estimates are shown in parentheses. 

The top spectrum displayed in Figure [5] was obtained 2007 July 5 at WIYN. In this 
spectrum, the secondary's spectral features (redshifted) are consistent with an 09III spectral 
type. Line ratios were estimated by a combination of equivalent width and line depth 
comparisons. The temperature-sensitive ratio, He I/He II A4026:He II A4200 A (~2:1), in 
addition to the presence of He I A4009 A, the ratio of He I A4144:He II A4200 A (1:2), 
and the ratio of He I A4471:Hy &: B .S (2:3) indicates the temperature class is most likely 
an 09 (IWalborn &: Fitzpatrickl Il990l ). The luminosity indi cative ratio He I A4026:Si IV 
A4089 A is ~2:1 and corresponds to a luminosity class of III (IWalborn fc Fitzpatrickl 119901 ) . 
The primary's features (blueshifted) are consistent with an 06V-III spectral type. The 
primary's spectrum contains little or no evidence of Si IV AA4089, 4116 A, Mg II A4481 A, or 
C III A4070 A indicating a temperature class 06 or earlier. The temperature-sensitive ratio, 
He I/He II A4026:He II A4200 A is ~1:1 and also indicates a temperature class of 06. It is 
difficult to distinguish He I A4009 A from the continuum noise, but there appears to be little 
or no evidence of He I A4144 A, and He I A4471 A is comparatively weak. Because Si IV 
absorption is weak even in an 06V-III and we did not include He II A4686 A or A4542 A 
in our spectral cover age, the luminosity c lass can only be constrained to lie between III and 
V (interpolated from IMartins et al.l 120051). The bottom spectrum in Figure is a composite 
of two spectra from the IWalborn & Fitzpatrickl Jl99ph digital atlas, HD37043 (09III) and 
HD93130 (06III). Intrinsic He I, He II, C III, Si IV, Mg II, and hydrogen absorption are 
labeled. We found the best agreement with Schulte 3 when HD37043 and HD93130 were 
combined with an equal component luminos ity ratio, indicating the primary may in fact have 
a luminosity class of IV (interpolated from IMartins et al.ll2005l ). 

There is a discrepancy between the theoretical mass ratio (q = 0.64 ± 0.02 for an 
06III & 09III or q = 0.72 ± 0.02 for an 06V & 09III) obtained from IMartins et al.l J2005h 
and the mass ratio of q = 0.44 ± 0.08 found in this study. We attribut e this discrepancy to 
either the uncertainty in mid-to-early O star masses (IMassey et al.ll2005l ). or the possibility of 
significant m ass transfer between the components based on th e hydrogen and he lium emission 
line profiles (IMiller et al.ll2007l ; iFalceta-Goncalvez et al.ll2006l ; ISana et al.ll200ll ) . It should be 
noted that an uncertainty in the spectral type of each component could also partially explain 
the discrepancy. Fo r instance, an O 611 1 primary and 09.5V secondary results in a mass ratio 
of q = 0.45 ± 0.02 (IMartins et al.l I2005T ) . However, this combination produces a luminosity 
ratio in disagreement with the implied luminosity ratio of the double-lined spectra and the 
light curve of Schulte 3. An 05III primary and 09.5111 secondary would also replicate the 
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computed mass ratio within the error (q = 0.51 ± 0.01; iMartins et al.ll2005l ). but also suffers 
from the same luminosity ratio disagreement. A likely scenario that explains the mass ratio 
discrepancy may in yolve a combination of two or more factors. Further analysis of the light 
curve for Schulte 3 (IKinemuchi et al.ll2008l in prep) and additional spectroscopic coverage is 
necessary however, to verify the presence of a mass transfer mechanism and determine what 
type of eclipsing system this is. 



5. Period Estimates of three Multi-lined Systems 

The 2007 August /September dataset was not sufficient to fit solutions to three double- 
lined systems but did allow for stricter constraints on the orbital elements. We applied the 
same method used for Schulte 3 to measure radial velocities for the components of these 
systems. Figures [6], HJ & M provide the time series for Ha A6563 A and He I A5876 A in 
velocity space for MT252, MT720, and MT771 (2007 August 28 through 2007 September 
4). Each system was introduced in Paper I with approximate periods and secondary mass 
estimates. 



5.1. The SB2 MT252 

For the first double-lined system, MT252, we obtained 23 observations from 1999-2007. 
MT252 is classified as a B1.5III in Paper I. We revise the classification slightly in this study 
to B2III. The primary shows weak C II/O II A4070 A, little or no Si IV A4089 A, and a 
ratio of Si III A4552:Si II AA4128, 4430 A of over 2:1, indicating a B2 temperature class. The 
luminosity indicator, Si III 4552:He I A4387 A is nearly 2:1 and points to a luminosity class 
of III. The secondary shows an absence of Si IV A4089 A and C II/O II A4070 A, making 
it no earlier than Bl. There is also no evidence of C II A4267 A or Si II AA4128, 4130 A, 
indicating it is likely no later than a B2. A He I A4121:He I A4144 A ratio of nearly 1:1 and 
the absence of N III A4634 A and O II AA4640, 4650 A absorption places the luminosity at 
V. We therefore conclude that the secondary is most likely a B1V. 

Of the 23 observations obtained, only six are of high enough S/N to estimate He I line 
velocities. These velocities listed in Table S] are relative to the heliocentric frame of rest. The 
largest line separation for MT252 is seen on 2007 August 28 (Figure [6]), and a near reversal of 
spectral features is seen within the ~8 days covered, indicating a likely period of 18-19 days. 
The maximum line separation in He I A5876 A is ~192 km s _1 , and the V r semi-amplitudes 
are K 2 = 107 ± 13 km s -1 and K\ = 85 ± 13 km s _1 relative to a systemic velocity of 
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7 = —18 ±10 km s 1 estimated from the 2007 August 30 and September 1 spectra (where the 
errors are approximated as they were with Schulte 3) indicating a mass r atio of q ~ 0.8 ±0.2, 



in ag reement with the theoretical mass ratio of q = 0.9 (interpolated from lDrilling &: Landolt 



20001 ). Based on a partial V r vs. time curve of this system, the eccentricity appears low. 



Therefore, for an assumed circular orbit, semi-amplitudes of K\ = 85 ± 13 km s 1 and 
K 2 = 107 ± 13 km s" 1 on 2007 August 28, and a period of P ~ 18.5 ± 0.5 days, the 
calculated masses are Misin 3 ~ 8 ± 1 M Q and M 2 sin 3 i ~ 6 ± 1 M . These values also imply 
a primary and secondary semi-major axes of aisin % ~ 22 ± 3 R Q and e^sin % ~ 27 ± 3 R Q . It 
should be noted that these values reflect lower limits in the absence of a full orbital solution. 



5.2. The SB2 MT720 

MT720 is composed of two components of apparently equal luminosity based on the 
similar depth and ratios of hydrogen and He I lines, as illustrated in a small portion of the 
2001 August 24 and 2001 September 9 observations in Figure [3 The top spectrum (August 
24) shows a nearly complete blend while the bottom spectrum (September 9) shows the 
deblended state (with perhaps a slight asymmetry in the red wing of the Balmer lines). The 
blended spectral features seen on 2007 September 1 & 3 (third and fifth spectra down) in 
Figure [8] indicate the systemic radial velocity, 7, is near zero. While the quality of the data 
makes it difficult to estimate an exact period, September 2 & 4 (fourth and sixth spectra 
down) indicate the period is ~ 5 days or less. The low S/N of the data makes it unclear 
whether the spectra show only two components. Therefore we applied both a two and three 
component model to both the hydrogen and He I lines in the attempt to deblend the spectral 
features. We obtained mixed results with the two-Gaussian fit (fixed widths), which provided 
velocity ratios varying from 0.4-0.9. Three Gaussians however, provided only slightly better 
fits to the data and again yielded varying velocity ratios. Therefore, we adopted the simpler 
two component model. The absence of He II in any of the spectra indicates that neither 
of the two components are earlier than about B0. The strength of He I coupled with the 
appearance of multi-component C III A4070 A and weak or absent Mg II A4481 A indicates 
that both are probably B2 or earlier. We therefore conclude the components lie between 
B0-B2 with uncertain luminosity classes. 



5.3. The SB2 MT771 



We obtained 36 observations (most low S/N) of the double system MT771. The primary 
is an 07V as indicated by the He II A4200:He I A4144 A and He II A4200:He I A4026 A ratios 
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of nearly 6:1 and 1:1 respectively. The ratio of He II A4542:He I A4387 A of just over 1:1 gives 
the luminosity class. The secondary is a probable 09V as best indicated by the temperature 
sensitive ratio of He II A4200:He I A4121 A (~1:1) and luminosity-sensitive ratio of Si IV 
A4116:He I A4121 A (just under 1:1). 

Figure [9] shows blended features for Ha A6563 A and He I A5876 A on days 2, 5, and 
8, suggesting periods of 1.5, 3, or 6 days. The similar He line profiles on days 1 and 3 
suggest that one full period has elapsed, leaving the ~1.5 day period as the only viable 
option. The August 30 (third spectrum down) spectrum shows a V r separation of 277 km s" 1 
( — 153 ± 14 km s _1 and 124 ± 13 km s _1 relative to the systemic velocity of 7 = —36 ± 
10.0 km s _1 estimated from the 2007 September 1 spectrum) and implies a mass ratio of 
q = 0.8±0.1, consistent with a th eoretical mass ratio of q = 0.68±0.02 for an 07V (25.9±0.6 
M ) and 09V (17.6 ± 0.5 M ) (IMartins et al.ll2005l ). Assuming a circular orbit, period of 

1.5 days, and minimum semi-amplitudes of K\ = 124 ± 13 km s _1 and K 2 = 153 ± 14 km s _1 
for the primary and secondary respectively, the minimum semi-major axes are aisin i ~ 

2.6 ± 0.3 R and a 2 sin i ~ 3.2 ± 0.3 R , and the masses are Misin 3 i ~ 1.8 ± 0.2 M and 
M2sin 3 z ~ 1.5 ± 0.2 M . As with MT252, these values reflect lower limits in the absence of 
a full orbital solution. No evidence for interacting winds or mass t ransfer is seen. Ther efore, 
using the theoretical radii of R\ = 9.3 ±0.1 and R2 = 7.6 ±0.1 (IMartins et al.ll2005l ). this 
system would need to have an inclination of i < 20° to avoid contact. This is consistent with 
the low inclination necessary (i ~ 24) to obtain the theoretical masses for an 07V and 09V. 



6. Summary of Survey Results to Date 

We have discovered the presence of six binary systems in the Cyg OB2 association 
as part of an ongoing survey to determine the distribution of binary orbital parameters. 
We presented the orbital solutions to three new systems (MT059, MT258, and Schulte 3) 
and constrained the orbital parameters to the three additional double-lined binaries (MT252, 
MT720 & MT771). Mass ratios for the six systems in this study range from q > 0.18±0.01 for 
MT258 to q ~ 0.8 ± 0.2 for MT252, and have periods ranging from ~1.5-19 days. Systemic 
velocities deviate by ~10-20 km s _1 from the cluster mean of —10.3 km s" 1 obtained in 
Paper I, but they are still well within the Association's radial velocity dispersion. These 
six systems more than double the number of known binaries in this cluster and bring the 
total number of uncovered OB binary/multiple systems in Cyg OB2 to 11. The locations of 
all 11 systems are shown in Figure [IOJ No evidence of grouping is apparent. Table [5] lists 
several key parameters for each system, including the star designation, photometric and/or 
spectroscopic binary type (Type), spectral classifications (S.C.), period (P), mass ratio when 
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available (q), and literature references (Ref). Periods for all 11 systems are relatively short 
(< 22 days) and the corresponding eccentricities are relatively small (0.03-0.11). With the 
exception of MT421, all systems have early companions. However, the apparent secondary 
mass distribution may depend on the sensitivity of the associated surveys. 

Although our current sample encompasses only a small fractio n of systems in Cyg OB2, 



we can test the limited mass ratio statistics against the hypothesis of lPinsonneault fc Stanek 



( 120061 ). that 45% of massive stars exist in "twin" systems with mass ratios of q > 0.95. Of 
the seven systems with mass ratio measurements in Table EJ we find at most one system that 
may be considered a possible "twin" (MT252 with q = 0.8 ± 0.2). If this hypothesis holds 
true, binomial statistics indicate the present distribution is likely only 8.7% of the time. On 
the other hand, if we consi der q < 0.95 for MT252, the p resent distribution is likely only 
1.5% of the time. Under the iPinsonneault fc Stanek I (120061 ) hypothesis, we would expect to 
observe a distribution of two to four "twins" in a sample of seven objects (21.4%, 29.2%, and 
23.9% respectively). Though no definite conclusions can yet be drawn about the agreement 
of these findings with IPinsonneault fc Stanek I (120061 ) owing to the limited number of systems 
available, at least in the case of Cyg OB2, a "twin" heavy binary distribution does not seem 
likely. 
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vatories for granting us observing time and making this project possible. This paper has 
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periences for Undergraduates (REU) program grant AST-0353760 and through grant AST- 
0307778, and the support of the Wyoming NASA Space Grant Consortium. The work of 
M. V. McSwain is supported by the National Science Foundation, under grant AST-0401460 
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Table 1. Summary of Observing Runs 



Date 


Observatory 


Spectral 


Grating 


Mean Spectral 


Date 




/Instrument 


Coverage 




Resolution 


Coverage 






(A) 


(1 mm -1 ) 


(A) 


(HJD) 


1999 Jul 4-5 


Keck/HIRES 


3890-6270 in 35 orders 


31.6 


0.1 


2,451,363-2,451,364 


1999 Jul 21-23 


Lick /Hamilton 


3650-7675 in 81 orders 


31.6 


0.1 


2,451,381-2,451,383 


1999 Aug 21-23 


Lick /Hamilton 


3650-7675 in 81 orders 


31.6 


0.1 


2,451,411-2,451,413 


1999 Oct 14-15 


Keck/HIRES 


3700-5250 in 29 orders 


31.6 


0.1 


2,451,466-2,451,467 


2000 Jul 10-11 


Lick/Hamilton 


3650-7675 in 81 orders 


31.6 


0.1 


2,451,736-2,451,737 


2000 Sep 18-19 


Keck/HIRES 


3700-5250 in 29 orders 


31.6 


0.1 


2,451,805-2,451,806 


2001 Aug 24 


WIYN/Hydra 


3800-4490 in order 2 


1200 


0.9 


2,452,146 


2001 Sep 8-9 


WIYN/Hydra 


3800-4490 in order 2 


1200 


0.9 


2,452,161-2,452,162 


2004 Nov 28-30 


WIYN/Hydra 


3800-4490 in order 2 


1200 


0.9 


2,453,338-2,453,340 


2005 Jul 18-21 


WIRO/WIRO-Spcc 


3800-4490 in order 1 


2400 


2.5 


2,453,570-2,453,573 


2005 Jul 18-20,22 


WIRO/WIRO-Spcc 


3800-4490 in order 1 


2400 


2.5 


2,453,632-2,453,635 


2005 Oct 13 


WIRO/Longslit 


4050-6050 in order 2 


600 


2.5 


2,453,657 


2006 Jun 16-20 


WIRO/Longslit 


3900-5900 in order 2 


600 


2.5 


2,453,903-2,453,907 


2006 Jul 15-16,20 


WIRO/Longslit 


3900-5900 in order 2 


600 


2.5 


2,453,932-2,453,935 


2006 Sep 8-11 


WIYN/Hydra 


3800-4490 in order 2 


1200 


0.9 


2,453,987-2,453,990 


2006 Oct 7 


WIRO/Longslit 


3900-5900 in order 2 


600 


2.5 


2,454,016 


2007 Jun 28-30 


WIRO/Longslit 


3900-5900 in order 2 


600 


2.5 


2,454,280-2,454,282 


2007 Jul 4-6 


WIYN/Hydra 


3820-4510 in order 2 


1200 


0.9 


2,454,285-2,454,287 


2007 Aug 28-Sep 4 


WIRO/Longslit 


5550-6850 in order 1 


1800 


1.5 


2,454,341-2,454,348 
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Table 2. Orbital Elements 



Element 


MT059 


MT258 


Schulte 3 


MT252 


MT771 


P (Days) 


4.8527 (0.0002) 


14.660 (0.002) 


4.7464 (0.0002) 


18-19 


1.5: 


e 


0.11 (0.04) 


0.03 (0.05) 


0.070 (0.009) 






u (deg) 


212 (17) 


68 (71) 


5.5 (0.7) 






7 (km s _1 ) 


-24.2 (1.7) 


-20.8 (1.2) 


-26.4 (1.7) 


-18 (10) 


-36 (10) 


T (HJD-2,400,000) 


53916.72 (0.22) 


53922.03 (2.87) 


53987.81 (0.01) 


54341.9 (1.0) 


54343.8 (0.5) 


K x (km s- 1 ) 


72.3 (2.3) 


41.2 (1.7) 


113.2 (14.5) 


85 (13) 


124 (13) 


if 2 (km s- 1 ) 






256.7 (2.4) 


107 (13) 


153 (14) 


M x (M Q ) 


21.4 (0.6) a 


21.4 (0.6) a 


> 17.2 (0.3) 


> 8± 1 


> 1.8 (0.2) 


M 2 (M ) 


5.1 (0.3)-13.8 (0.6) 


4.1 (0.2)-17.5 (1.0) 


> 7.6 (1.4) 


> 6±1 


> 1.5 (0.2) 


/(m)i (M ) 


0.187 (0.018) 


0.106 (0.013) 


0.709 (0.273) 






/(m) 2 (M ) 






8.289 (0.237) 






S. C.i 


08V 


08V 


06IV: 


B1.5III 


07V 


S. C.2 


B 


B 


09III 


B1V 


09V 


aisin i (Rq) 


6.9 (0.2) 


11.9 (0.5) 


7.4 (0.9) 


22 (3) 


2.6 (0.3) 


a2Sin I (Rq) 






16.7 (0.2) 


27 (3) 


3.2 (0.3) 


rmsx (km s - 1 ) 


11.0 


6.7 


38.7 






rms2 (km s - 1 ) 






9.0 







Note. — Calculated errors are located in parentheses. 
a Theoretical masses adopted from lMartins et all (|2005h . 
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Table 3. Ephemerides for MT059 & MT258 



Star 


Date 
(HJD-2,400,000) 




V r 

(km s- 1 ) 


l(j err 

(km s- 1 ) 


O-C 
(km s- 1 ) 


MT059... 


51381.50 


0.562 


10.0 


25.6 


-5.8 




51411.50 


0.744 


-89.1 


31.6 


-37.2 




51413.50 


0.156 


-4.0 


15.4 


15.4 




51467.91 


0.369 


51.8 


13.4 


11.2 




51736.50 


0.717 


-32.4 


32.8 


8.2 




51737.50 


0.924 


-98.0 


14.6 


5.2 




51805.75 


0.988 


-106.7 


6.0 


-10.5 




52146.70 


0.248 


15.7 


4.1 


-2.7 




52162.66 


0.538 


26.9 


4.4 


4.7 




53338.62 


0.871 


-90.3 


4.7 


6.1 




53570.50 


0.654 


-26.3 


5.2 


-11.3 




53571.50 


0.860 


-103.3 


5.2 


-9.4 




53573.50 


0.272 


21.3 


5.1 


-4.3 




53633.50 


0.637 


-15.6 


4.9 


-7.1 




53634.50 


0.843 


-88.6 


5.0 


0.4 




53636.50 


0.255 


26.4 


5.6 


5.6 




53657.50 


0.582 


34.7 


4.3 


24.7 




53903.72 


0.322 


28.3 


15.6 


-7.6 




53903.73 


0.324 


24.3 


11.1 


-11.9 




53903.74 


0.325 


32.6 


11.0 


-3.7 




53903.84 


0.346 


37.0 


11.5 


-1.8 




53903.84 


0.347 


39.1 


12.0 


0.0 




53904.70 


0.523 


5.7 


11.6 


-20.0 




53904.71 


0.525 


35.0 


11.5 


9.6 




53904.86 


0.558 


7.0 


11.1 


-10.1 




53904.89 


0.563 


22.9 


11.0 


7.1 




53905.74 


0.739 


-82.5 


13.8 


-33.0 




53905.82 


0.755 


-53.0 


11.0 


3.4 




53905.84 


0.758 


-48.6 


10.9 


9.1 




53906.72 


0.939 


-99.7 


12.1 


3.4 




53906.72 


0.941 


-100.0 


11.4 


3.1 




53906.87 


0.972 


-112.8 


11.7 


-13.1 




53907.77 


0.156 


-22.4 


11.1 


-2.7 




53907.86 


0.175 


-10.3 


11.0 


0.0 




53932.73 


0.301 


42.1 


11.0 


9.9 




53932.83 


0.321 


77.8 


13.3 


42.0 




53932.84 


0.322 


67.8 


12.4 


31.8 




53935.77 


0.927 


-91.2 


13.2 


12.1 




53935.85 


0.943 


-71.5 


12.0 


31.4 




53987.70 


0.628 


-1.2 


8.3 


3.8 




53988.74 


0.843 


-79.0 


6.3 


10.0 




53989.66 


0.031 


-82.7 


4.2 


-0.9 




53990.85 


0.277 


30.2 


5.1 


3.2 




54285.93 


0.084 


-48.9 


6.4 


8.0 


MT258... 


51467.91 


0.262 


-11.4 


7.7 


1.3 
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Table 3 — Continued 



Star Date 

(HJD-2,400,000) 


<P 


V r 

(km s- 1 ) 


lcr err 

(km s- 1 ) 


O-C 
(km s- 1 ) 


51805.74 


0.223 


1.4 


4.8 


3.5 


52146.70 


0.386 


21.7 


3.9 


6.0 


52161.81 


0.368 


12.2 


3.8 


1.3 


52162.66 


0.424 


-1.7 


3.8 


0.1 


53338.62 


0.867 


-44.4 


3.6 


9.7 


53340.59 


0.995 


-57.4 


4.0 


3.0 


53571.50 


0.005 


-36.3 


4.4 


-7.1 


53573.50 


0.135 


-57.2 


4.5 


-0.1 


53632.50 


0.970 


-62.0 


6.3 


-2.4 


53633.50 


0.035 


-72.8 


4.7 


-11.5 


53903.76 


0.603 


12.5 


11.1 


-4.8 


53903.77 


0.604 


13.3 


11.3 


-4.1 


53904.76 


0.668 


11.9 


11.0 


-9.0 


539U4. / / 


0.668 


10.3 


11.6 


-10.5 


53906.75 


0.797 


-4.1 


11.0 


-9.0 


53906.76 


0.798 


1.7 


11.2 


-3.1 


53906.76 


0.798 


-5.6 


11.1 


-10.3 


53907.81 


0.867 


-17.7 


10.7 


-4.5 


53907.82 


0.867 


-12.6 


10.8 


0.6 


53932.75 


0.488 


11.3 


10.9 


-3.3 


53932.86 


0.495 


34.9 


11.1 


19.2 


53935.81 


0.686 


5.3 


12.6 


-3.3 


53935.86 


0.690 


19.6 


11.5 


11.6 


53987.70 


0.060 


-38.5 


4.6 


1.5 


53988.74 


0.127 


-17.0 


4.4 


6.5 


53989.66 


0.187 


-13.2 


3.7 


-4.7 


53989.78 


0.195 


-11.3 


4.2 


-4.7 


53990.85 


0.264 


8.2 


4.3 


-0.6 


54285.93 


0.446 


19.3 


4.0 


-1.6 


54286.66 


0.493 


21.2 


5.7 


2.7 


54342.88 


0.708 


10.9 


13.6 


-0.1 


54343.76 


0.768 


23.9 


14.0 


5.1 


54345.82 


0.908 


22.0 


13.8 


7.7 


54347.81 


0.044 


-6.7 


13.6 


10.0 
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Table 4. Ephemerides for Schulte 3, MT252, k MT771 



Star 


Date 




V r i 


Oi 


-Ci 


Vr2 


2 - C*2 




(HJD-2,400,000) 




(km s- 1 ) 


(km 




(km s- 1 ) 


(km s- 1 ) 


Schulte 3 


53987.69 


0.974 


-111 (20) 




35.4 


250 (8) 


2.9 




53988.74 


0.195 


12 (20) 




63.7 


13 (8) 


-7.7 




53989.65 


0.387 


78 (20) 




17.0 


-224 (8) 


6.1 




53989.77 


0.412 


46 (20) 




-22.6 


-244 (8) 


2.2 




53990.85 


0.639 


90 (20) 




47.9 


-179 (8) 


-5.1 




54285.93 


0.809 


-123 (20) 




-56.1 


91 (8) 


15.8 




54341.80 


0.580 


112 (20) 




46.7 


-229 (8) 


0.3 




54342.72 


0.773 


-10 (20) 




31.6 


12 (8) 


-6.6 




54343.83 


0.007 


-139 (20) 




7.7 


243 (8) 


-2.2 




54344.76 


0.203 


-20 (20) 




25.5 


13 (8) 


6.5 




54345.78 


0.418 


73 (20) 




2.7 


-257 (8) 


-7.4 




54346.79 


0.631 


69 (20) 




23.1 


-179 (8) 


3.9 




54347.75 


0.833 


-41 (20) 




42.9 


103 (8) 


-10.3 




54348.84 


0.063 


-125 (20) 




8.0 


211 (8) 


0.9 


MT252 


54341.85 
54342.84 
54343.79 
54345.85 
54346.88 
54348.78 




67 (8) 
1: 
-9: 
-27: 
-63: 
-64: 






-125 (8) 
-83: 

41: 




MT771 


54341.81 
54342.81 
54343.82 
54345.83 
54347.82 
54348.77 




63 (10) 

88 (8) 
-36 (10) 






-165 (10) 
-11: 
-189 (10) 
-36 (10) 
89: 
-26: 
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Table 5. OB Binaries in Cyg OB2 



Star 


Type 


S.C. 


P 

(days) 


q 


Ref. 


MT059 


SB1 


08V & B 


4.8527 (0.0002) 


0.24 (0.02)-0.64 (0.03) 


1 


MT252 


SB2 


B2III & B1V 


18-19 


0.8 (0.02) 


1 


MT258 


SB1 


08V & B 


14.660 (0.002) 


0.19 (0.01)-0.84 (0.05) 


1 


MT421 


EA 


09V & B9V-A0V 


4.161 




2 


MT429 


EA 


BOV & ?? 


2.9788 




2 


MT696 


EW/KE 


09.5V & early B 


1.46 




3 


MT720 


SB2 


early B & early B 


< 5 




1 


MT771 


SB2 


07V & 09V 


1.5: 


0.8 (0.1) 


1 


Schultc 3 


SB2/EA: 


06IV: & 09III 


4.7464 (0.0002) 


0.44 (0.08) 


1,4 


Schulte 5 


EB 


07Ianfp & Ofpe/WN9 (& BOV?) 


6.6 


0.28 (0.02) 


5,6,7,8,9,10 


Schulte 8a 


SB2 


05.51 & 06: 


21.908 


0.86 (0.04) 


11,12 



Note. — Photometric types EW/KE, EA, and EB stand for Contact system of the W UMa type (ellipsoidal; P < 1 
day), Algol type (near spherical), and j3 Lyr type (ellipsoidal; P > 1 day) respectively. 
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Fig. 1. — K curve and orbital solution to the 08V star, MT059. The symbols correspond 
to the different observatories, where the stars are observations taken with the Hamilton 
Spectrograph (Lick), the diamonds are observations taken with HIRES (Keck), the triangles 
are observations taken with Hydra (WIYN), and the squares are observations taken with 
WIROspec or the Longslit Spectrograph (WIRO). 
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Fig. 2. — V r curve and orbital solution to the 08V star, MT258 using the same format as 
Figure [Q 
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Fig. 3.— Ha A6562.80 A (left) and He I A5875.75 A (right) in velocity space for Schulte 3 
between 2007 August 28 and 2007 September 4. Dates are listed as HJD - 2,400,000 and 
apply to both left and right panels. 
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Fig. 4. — V r curve and orbital solution for Schulte 3 using 14 of the 18 observations. The 
filled points correspond to the primary (06IV:) and the unfilled points correspond to the 
secondary (09III). The symbols represent the same observation locations as Figure [H 




Fig. 5. — A spectrum of Schult e 3 obtained on 2007 July 5 wit h WIYN (top), and a composite 
of two shifted spectra from the lWalborn fc Fitzpatrickl (1l990l ) digital atlas (bottom) that best 
match the spectral types of the primary and secondary of Schulte 3. The primary is shifted 
-110 km s _1 and the secondary is shifted 257 km s -1 . Absorption features are labeled at 
their respective intrinsic values. 
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Fig. 6. — Ha A6562.80 A and He I A5875.75 A absorption in velocity space for MT252 using 
same format as Figure [31 



-27- 




4300 4350 4400 4450 4500 

Wavelength (A) 

Fig. 7. — A portion of the 2001 August 24 (top) and 2001 September 9 (bottom) spectra for 
MT720, showing the nearly blended and unblended states respectively. 
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Fig. 8. — Ha A6562.80 A and He I A5875.75 A absorption in velocity space for MT720 using 
same format as Figure [31 
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Fig. 9. — H« A6562.80 A and He I A5875.75 A absorption in velocity space for MT771 using 
same format as Figure [31 
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Fig. 10. — A map of OB stars from this survey. Circles indicate the location of known OB 
binaries within this region of Cyg OB2. 



